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Abstract. In a 2D rectangular microwave cavity dressed with one point-like 
scattcrer, a semiclassical approach is used to analyze the spectrum in terms 
of periodic orbits and diffractive orbits. We show, both numerically and 
experimentally, how the latter can be accounted for in the so-called length 
spectrum which is retrieved from 2-point correlations of a finite range frequency 
spectrum. Beyond its fundamental interest, this first experimental evidence of 
the role played by diffractive orbits in the spectrum of an actual cavity, can be 
the first step towards a novel technique to detect and track small defects in wave 
cavities. 



PACS numbers: 05.45.Mt, 03.65.Sq 
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1. Introduction 

It has been widely recognized that the semiclassical approach of spectral statistics 
based on periodic orbits (POs) has met a definite success in chaotic systems (see 
e.g. [P and references therein). In systems where one or more point-like scatterers 
are added, a similar approach is possible, calling for both POs and the so-called 
diffractive orbits (DOs) While the proliferation with length of diffractive orbits 
may be reminiscent of the proliferation of POs in chaotic systems, if the unperturbed 
system is regular, the statistics is clearly not predicted by Random Matrix Theory 
even if some level repulsion appears due to the presence of scatterers. Indeed, at large 
spacing, the spacing distribution decays exponentially like in a regular system 02]. 
In fully chaotic systems, on the contrary, adding scatterers has been shown to be of 
practically no effect on the spectral statistics 

In the present paper, we will show both numerically and experimentally, how the 
presence of short DOs can be accounted for in the so-called length spectrum which 
one can retrieve from two-point correlations of a finite range frequency spectrum. 
Our experiments are performed in microwave cavities at room temperature and the 
limited range of the available spectrum is a practical issue. Nevertheless, in the case of 
a 2D rectangular (regular) cavity with one point-scatterer (actually a small scatterer 
in the experiments), the lengths associated to the shortest DOs are easily identified 
in our experiments and can be used to locate the scatterer and possibly track its 
displacements. 

2. Semiclassical expansion and spectral density 

We now briefly recall the semiclassical approach of spectral correlations based on 
a decomposition of the spectral density on POs. It relics on the semiclassical 
representation of the Green's function G(f,f'), written as a sum over classical 
trajectories connecting points r and r ' 

G(r, r') = J2 Atr{r, 7) exp (iS tr (r, r ')) (1) 

cl.tr. 

By evaluating the trace of G{f, r) through a stationary phase approximation, one 
obtains a semiclassical trace formula for the modal density p pQ. The latter is thus 
written as an average smooth density plus an oscillatory part: 

p(k) = p(k) + p° sc (k) = p(k) + £ + c.c. (2) 

i 

Here, k is the wavenumber, £j is the length of the jth PO and Aj is its complex 
amplitude accounting for its stability and possibly depending on k. 

In diffractive systems with point-like singularities, classical trajectories that hit 
those singularities can be continued in any direction. These can nonetheless be tackled 
with in the wave description by introducing an isotropic diffraction coefficient T>, which 
fixes the scattering amplitude at each scatterer. In previous works (see [01 El), this 
diffraction constant has been calculated (with the free Green's function in [Tj) to yield: 
_ 2tt 

V = -ln(fea/2) - 7 + i7r/2 (3) 
where 7 is the Euler constant and a is a characteristic length which may be interpreted 
as the non- vanishing radius of an s-wave scattering disc [H], the above expression being 
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precisely the limiting scattering amplitude for ka <C 1. The semiclassical expansion of 
the Green's function in the presence of a point scatterer located at s therefore reads: 

G(r, r') = G Q (r, r') + G (r, s) V G(s, r') (4) 

where Go is the unperturbed Green's function. Hence formula J5J for the modal 
density still holds provided that nonclassical contributions due to DOs arc included, 
consistently with the geometrical theory of diffraction, yielding: p osc = p°^ + p d s c . 
In a rectangular domain of area A with a single point scatterer, contributions from 
periodic orbits and diffractive orbits respectively read [S]: 

A 1 00 k 

P°p S o C (k) = -J2 ,o n C0 < rM v° - rn po* - (5) 

and 

P° d o(k) = Y! — -JL= cos(k£ do - n do n - 3tt/4) (6) 

where denotes a sum over primitive periodic (diffractive) orbits of length £ po (£do) 
and number of bounces n po (rido), T is the number of repetitions. In formula © 
only leading order one-scattering events are included, repetitions or concatenations of 
primitive orbits of order v being of order 0{k~ y l 2 ) [S]. 

Note that the contribution of DOs is subdominant with respect to POs due to 
different fc-dependences. Nevertheless, the relevance of both POs and DOs is clearly 
illustrated through a weighted length spectrum |5] which is obtained by calculating the 
so-called form factor K(L), i.e. the Fourier transform of the spectral autocorrelation 
C(k) oip osc {k): 

C(K) = (p° sc (k + ^)p° sc (k-^) k (7) 

In practice, the local average over k in JJJ) can be written 

(f(k)) k = J dk' f{k')W a {k' -k) (8) 

where the weighing function W a is zero-centered and of typical width a. In the 
following, we will use either Gaussian or Hanning weighing functions. 

In the so-called diagonal approximation , the following expression of the form 
factor is obtained |l()j : 

KWoc^AjfSiL-ej) (9) 

3 

Ideally, the length spectrum appears as a series of delta peaks located at particular 
orbits lengths ii with real positive amplitudes (see Figure QJ. 

For practical reasons, both in numerics and in analyzing our experiments, we will 
rather use the cumulated density N(k) = dk' p(k'), which is a staircase function 
increasing of one unit at each eigenwavenumber. This integrated quantity enables one 
to evaluate the form factor more readily. 
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Figure 1. Length spectrum computed in a rectangular cavity with a single point 
scattcrcr. Approximately 11000 resonances have been used. Green sticks indicate 
the lengths of the POs listed in Table Q] 



3. Numerical length spectra 

To validate the possiblity of identifying the shortest DOs in the short-term non- 
universal part of the length spectrum, we first investigate a numerical model of a 
rectangular cavity with a point scatterer (hereafter called the dressed cavity). To 
numerically obtain all the eigenfrequencies in a given frequency range for the dressed 
cavity, we will closely follow the method first introduced by R. L. Weaver and D. 
Sornctte in reference This method provides a clear physical insight of the 

way the scatterer can be characterized by its transition matrix and scattering cross 
section. However, following a criticism made by Shigehara and Cheon in we wm 
correct a pathology of the aforementioned method by restoring the orthogonality of 
the perturbed cigenfunctions. 

In 2D billiards, in the presence of a point scatterer at position s, it is possible to 
write the total field "J at the position r, in the vicinity of s, as the superposition of 
incoming and outgoing waves: 



l -H^\k\f- s\) + Q - i|) H^\k\f- s\) 



(10) 



(1 2) 

where Hq ' arc the Hankcl functions of first and second kinds. The parameter t(k) 
defines the ratio of the strength of the incident field at the scatterer to the strength 
of the outgoing field in the vicinity of the scatterer. In terms of this parameter the 
scattering cross section (here a length) is given by |£| 2 /4fc. From the expression (|10fl 
for 'J, energy conservation implies | ^ | 2 =| h — ij | 2 yielding t = (a + i/4) — 1 , a being 
a real parameter. 

In the appendix, we show how a specific fc-dependence of a can be established 
to ensure the orthogonality of the perturbed cigenfunctions. We thus reconcile the 
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seemingly contradictory approaches of Refs [H] and ^21 • The resulting expression for 
t (equation ljA.15| l1 is the same as the diffraction constant in equation ©. 

Note that, apart from a logarithmic correction, the scattering cross section 
essentially scales as the wavelength, thus making the scattercr equally efficient at 
all frequencies, which is not the behavior of finite size scatterers as used in our 
experiments. However, in practice, we will always be in the limiting case ka <C 1. 
The resonances of the dressed cavity correspond to the poles of the transition matrix 
l)A.2|) . Practically, the ensuing characteristic equation is solved through a convergence 
accelerating procedure described in JT], by equating a(k) = — [ln(fca/2) +7] /(2n) 
to the function g(k): 

g(fe) = E^ fc8 _ (fc co) )8 -4 < u > 

Here, the ^„'s are the solution of the Helmholtz equation and the fc^ c) 's are the 
associated eigenwavenumbers in the empty cavity. The problem of calculating the 
cigenwavenumbers in the presence of the point scatterer is thus solved by finding the 
zeroes of g(k) — a(k). Note that g(k) — a(k) has singularities for each kn ] and all 
its zeroes, hereafter denoted k i: must lie between two consecutive eigenvalues of the 
empty cavity. This analytical approach allows us to deal with very large frequency 
ranges. 

In practice, we obtain the form factor and the length spectrum by calculating the 
Fourier transform of the difference N(k) — N(k), where the average behavior N(k) is 
given by Weyl's formula: 

N( k ) = ^k 2 -^k+\ (12) 

47T 47T 4 

Here A is the total area of the cavity and C its perimeter. 

An example of a length spectrum corresponding to a rectangular cavity with 
a single point scatterer is given in Figure ^ Here, the dimensions of the cavity in 
which we have performed our calculations are those of the actual cavity we use in 
our experiments : L x = 0.7562 m and L y = 0.4656 m, the scatterer is located at 
x = 46.5 cm and y = 20 cm (positions are measured from the upper left corner of the 
cavity sketched in Figures [5] and 0} , and the value of a = 3 mm corresponds to the 
radius of the actual cylindrical scatterer. A large number of resonances (approximately 
11000) have been used so that the length resolution is excellent. To illustrate that 
such a length spectrum still is dominated by the POs of the empty cavity, it is shown 
in Figure ^ with an amplitude scale such that the contributions of the DOs are too 
small to be seen. 

In a rectangular cavity of side lengths L x and L y , the POs are identified by 
two indices n and m indicating the number of steps on a rectangular lattice of basis 
(2L x ,2L y ) (see reference 53)- They form families of continuously deformable orbits 
with the same length £ = 2(n 2 L 2 + m 2 -^ 2 ) 1 ' 2 . For the families shown in Figure [21 
the most symmetric representative is displayed, as well as the result of continuous 
deformations (red or blue) yielding self-retracing orbits hitting two corners. Note 
that, in polygonal billiards, diffraction occurs only at vertices where the angle is not a 
submultiple of it [S]. POs arc easily identified on the length spectrum shown in Figure 
njwith the help of Tablc^Jwhich gives the correspondence between indices and lengths 
for all lengths shorter than 5 meters. At first sight, it could even seem that no other 
contribution is to be seen as if the DOs were absent from it. Somehow, it could even 




Figure 2. Periodic orbits representation with indices (a) (2,3) (b) (1,4). 
Continuous deformations of the orbit shown in thick black lines lead to the orbits 
shown in red or blue lines 



be expected since no long range correlations arc observed in the frequency spectrum 
thereby indicating that if the DOs should contribute, especially at short lengths, they 
should in a negligible way. 



Table 1. Indices and lengths (in meter) of the periodic orbits shown in Fig. 1 



(n, m) 


length 


(n, m) 


length 


(n, m) 


length 


(n, m) 


length 


(0,1) 


0.931 


(0,3) 


2.794 


(0,4) 


3.725 


(0,5) 


4.656 


(1,0) 


1.512 


(2,0) 


3.025 


(1,4) 


4.020 


(2,4) 


4.798 


(1,1) 


1.776 


(2,1) 


3.165 


(2,3) 


4.118 


(1,5) 


4.896 


(0,2) 


1.862 


(1,3) 


3.177 


(3,0) 


4.537 


(3,2) 


4.905 


(1,2) 


2.399 


(2,2) 


3.552 


(3,1) 


4.632 







This is indeed what one can observe by closely inspecting a typical length 
spectrum for lengths smaller or of the order of the size of the cavity in the presence 
of a single point scatterer. In Figure [31 the contributions of the DOs are displayed 
on the length range from to 1.6 m using an enlarged scale for the amplitude of the 
peaks. Sticks with different colors indicate the lengths of the DOs (blue) and POs 
(green) within this range. 

The shortest relevant lengths are those of the elementary self-retracing DOs 
starting from the scatterer and bouncing once on one of the sides of the rectangle 
as shown in Figure |I{a). In Figure 0[b), some of the shortest self-retracing DOs with 
one scattering event, two bounces on one side and one bounce in a corner are drawn. 
In Figure0{c), DOs with one scattering event and three bounces are drawn. All these 
DOs are listed with their corresponding lengths in Table [3 

In Figure El we have used different Gaussian normalized frequency windows with 
variance a = 2 GHz centered on 10, 15 and 20 GHz to illustrate the l//c-dependence 
of the squared amplitudes associated to the DOs (see formula ©). An excellent 
agreement with the prediction of formula © is observed and validates the semiclassical 
approach. We also have observed the fc-dependence expected for POs (see formula (0)), 
but we will not present it here. 

In this section we have used a model for a point-scatterer in a cavity to validate 
the possibility of revealing the presence of DOs in the short-term non-universal part 
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Figure 3. Zoom (xlO 4 ) of the length spectrum shown in Figure HI on the length 
range from to 1.6 m using an enlarged scale for the amplitude of the peaks : 
POs (green sticks), DOs (blue sticks) 



Table 2. Lengths (in meter) of the diffractive orbits shown in Figure 4; 1st, 2nd 
and 3rd kinds refer to DOs depicted in (a), (b) and (c) respectively. 





1st kind 




2nd kind 


3rd kind 


a 


0.400 e : 


0.789 


a : 1.454 


a : 1.099 


1) 


0.532 f : 


0.930 


b : 1.575 


b : 1.316 


c 


0.582 g : 


1.012 


c : 1.625 


c : 1.564 


d 


0.706 h : 


1.071 


d : 1.734 


d : 1.603 



of the length spectrum. 

4. Experimental length spectra 

In our experiments, the frequency spectra are determined from transmission signals 
measured in a 2D microwave cavity operated at frequencies ranging from 500 MHz to 
5 GHz. The rectangular cavity is composed of two copper plates sandwiching a copper 
rectangular frame of thickness 5 mm and of width 2 cm. The cavity may thus be viewed 
as the slice of a rectangular waveguide closed at both ends, with perimeter C = 2.446 m 
and area A = 0.3528 m 2 . The quality of copper is OFHC to reduce ohmic losses. Due 
to its height of 5mm (smaller than half the smallest wavelength used), this cavity 
only admits transverse magnetic two-dimensional modes of order 0. Through one of 
the copper plates, ten antennas are introduced with positions determined at random. 
Optimal coupling was obtained by fixing their penetration length inside the cavity 
at 2 mm. The details about transmission measurements are described in reference 
[Tl] . The scatterer is a small copper cylinder of radius 3 mm much smaller than the 




Figure 4. Three different kinds of one-diffraction event diffractive orbits 
associated to the position of the scatterer used numerically and experimentally: 
(a) direct self-retracing DOs (1st kind), (b) self-retracing DOs (2nd kind), (c) 
3-bounce DOs (3rd kind). 
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Length (m.) 

Figure 5. Scaling law in 1/k of the amplitudes of the peaks associated to the DOs. 
Gaussian normalized frequency windows with equal variances (2 GHz) centered 
on 10 GHz (green), 15 GHz (blue) and 20 GHz (black) were used. An excellent 
agreement with the prediction of formula J^J is observed. 



smaller wavelength of the order of 5 cm. For the results presented here, we used three 
different couples of antennas to be sure not to miss any resonance frequency in the 
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dressed cavity in the range mentioned above (nearly 300 resonances were measured 
for each position of the scattcrcr). The cumulated density number N(k) we could 
deduce from these measurements enabled us to verify that a small level repulsion can 
be observed as exemplified by the nearest-spacing distribution P(s) whose histogram 
is shown in Figure The nearest-spacings are as usually obtained by building the 
sequence of normalized spacings s, = N(ki+i) — N(ki), whose average is unity. A 
comparison is shown with the distribution associated to a semi-Poisson sequence as 
a guideline. In fact, Bogomolny et al. have shown that the statistics of the singular 
billiard as the one discussed in the previous section, is intermediate between Poisson 
(uncorrelated spectrum) and semi-Poisson (short range repulsion) 

0.25 

0.2 
0.15 

0.1 
0.05 


1 2 3 4 5 



Figure 6. Experimental histogram of P(s) in the dressed cavity; Semi-Poisson 
law: P(s) = 4se — 2s (blue curve). 

Once the frequencies of all the resonances in the range mentioned above have been 
determined with a sufficient precision (better than one tenth of the local spacing), 
we build an average cumulated density through a polynomial fitting of order two 
(consistent with Weyl's formula). Hence, an experimental length spectrum is evaluated 
by using the total frequency range at our disposal (here, the weighing function of 
formula (JSJ) is therefore simply a Hanning window) for the FFT. Figure shows such 
a length spectrum on a scale where the shortest POs can clearly be identified. Note 
however the poor resolution compared to the analytical length spectrum of Figures |5] 
or|SJ This is entirely due to the reduced frequency range we were compelled to use. 
Indeed, above 5 GHz, the ohmic losses are so important that the overlap of neighboring 
resonances prevents from properly extracting all the resonance frequencies with the 
required precision. If one wishes to use frequencies above this limit, the resolution of 
the length spectrum is not significantly improved whereas the risk of missing many 
levels is rapidly increasing, thereby making the appearance of spurious peaks in the 
length spectrum more probable. These arc current limitations when one requires high 
precision spectra in such microwave experiments at room temperature. Within this 
frequency range the resolution is sufficient to identify the shortest (and relevant) DOs, 
even when only one couple of antennas is used to extract the frequency spectrum. 

As in the analytical part, one can note the presence of peaks not centered at the 
lengths of POs. The important difference with the analytical or numerical approach is 
that the peaks related to the DOs have amplitudes only one order smaller than those 
of the POs. This can be accounted by the finite size of the actual scatterer leading 
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Figure 7. Experimental length spectrum computed from 300 actual resonances 
measured from 500 MHz to 5 GHz. The position of the center of the scatterer is 
the same as in the numerical results of the previous section, (a) length range from 
to 5 m, POs indicated by green sticks; (b) length range from to 1.6 m, POs 
indicated by green sticks, DOs by blue sticks. 



to a different fc-dependence of the amplitudes associated to the DOs. In addition, 
we have also checked that the displacements of the scatterer within the cavity can 
be tracked by properly monitoring the lengths of the peaks associated to the shortest 
DOs contributing to the experimental length spectrum. This is illustrated in Figure 
[H] where only one couple of antennas was used to extract the frequency spectrum for 
four different positions of the scatterer. This could constitute the basis of a novel 
non-destructive way of detecting the motion of a defect. 

5. Conclusion 

In the present paper, we have shown how a semiclassical approach of spectral statistics, 
usually based on periodic orbits in chaotic cavities, can be extended to analyze 
experiments in a 2D rectangular microwave cavity with a small scatterer by including 
diffractive orbits. More specifically, we have illustrated, through a model of cavity with 
a point-scatterer, how DOs can be clearly identified in the so-called length spectrum 
of the cavity. This length spectrum could be retrieved through the evaluation of the 
two-point spectral correlations deduced from the cumulated spectral density of around 
300 resonances measured in the frequency range [500 MHz, 5 GHz]. To our knowledge, 
this is the first experimental evidence of the unambiguous role played by diffractive 
orbits in the spectrum of an actual cavity. This is also a clear illustration of how 
a semi-classical analysis can be used to reveal non-universal spectral features of a 
complex cavity. Our findings show that the length spectrum can be used to locate the 
scatterer, or any small defect, and, eventually, to track a moving defect. 
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Figure 8. Experimental length spectrum for 4 different positions of the scattcrcr 
measured from the upper left corner of the rectangular cavity: (a) x = 28.5 cm, 
y = 13 cm, (b) x = 43.5 cm, y = 23 cm, (c) x = 45 cm, y = 23 cm, (d) 
x = 46, 5 cm, y = 23 cm. The corresponding DOs of the three kinds are depicted 
by blue sticks. POs are indicated by green sticks. 
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Appendix A. Orthogonality of the perturbed eigenfunctions 

We will use the following notations: Go is the Green function of the bare cavity, ^ n 
axe the eigenfunctions and E n the corresponding eigenenergics; G is the perturbed 
Green function of the dressed cavity and denotes the eigenfunctions with z n the 
corresponding eigenvalues. The two Green functions are related by: 

G (r, r"; z) = G (r, f; z) + G Q (r, s; z) t (s; z) G (s, f '; z) (A.l) 

where the transition matrix r is given by 

r(s;z)=[t- 1 -f(s;z)}' 1 (A.2) 
with t^ 1 = a + i/4 (see formula (12) of reference ^J) for some real a. 
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The resonances of the perturbed cavity correspond to the zeroes of [t 1 — / (s; z)] 
where 



f{s;z) = lim 



r — >s 

i 



Go + (VZ\\r-s\\) 



, lim 

4 r— >s 



G (f, s; z) - — ( 7 + In 



V5||r- 



One deduces the 4>„ through the residues of G at z„ 

lim [G (f, s; z) (z - z„)] = $„ (f) (s) 

z — >z n 

Thus, according to (|A.1|I (see 22); one has, 
$„ (f) = 7V„G (f,s;z n ) 

where 

^(z„-iS„)V 
One wishes to calculate / ^> m (f) ^> n (f) dr. Using 

*v (f) (s) 



Go (f, s;z) = ^2 



z-E v 



one gets, 



<t> m (r) (f) df 

(«) *m (5) 



= N m N n Y, 



(z m - (z„ - 



ifi^fr] df 



- N N T ^ 

where J ^ v (f) ^ (f) df = (5,,^ was used. If m = n, one finds: 



[<t>l(r)df=N*J2 7 1 

J Z \ Z n 



(z„ - E u y 



If n 7^ m, one calculates 

/ (s; z n ) - / (s; z m ) 



lim 

r — >s* 

E 



G (f, s; z„) - G (f s; z m ) - — In 

47T 



(2n - E v ) 
fern Zn) ^ ^ 



E 



in 



(z m -Ev) 4n~\z. 
*l(s) 1 



In 



(z„ - E v ) (z m - E v ) 47r \z. 



(A.3) 
(A.4) 
(A.5) 

(A.6) 
(A.7) 



(A.8) 



(A.9) 



(A.10) 



At the resonances of G, / (s) = t 1 assumed to be constant, whence / (s, z 
/ (s, z m ), thus implying, if one writes z n = k^ : 

1 



E 



(z„ - Ev) (z m - Ev) 2n(kl-k n ) 



(A.11) 
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If one allows a to depend on k, then one can restore the orthogonality between the 
eigenfunctions <& n and <& m with m ^ n. Indeed, it suffices to write: 

a(k n )-a(k m ) = -^-ln(^-) (A.13) 
2tt \ k m J 

or, equivalcntly, 

a(k) = — — In ( ) + const (A. 14) 

2n \ k J 

If one puts fco = 2/a and const = — 7/(27r), 7 being the Euler constant, then (see 

*= nrr^ ( A -!5) 

-ln(f)- 7 + if 

This corresponds to the asymptotic behavior of the diffraction coefficient of a 
cylindrical s-wave scattercr of radius a. Therefore, in the approach due to R. L. 
Weaver and D. Sornctte [TT], the relation (23) which defines the resonances should 
read: 

- s h (t) - fc (G » « »• => - s k (^) ' (A - 16) 

thus restoring the orthogonality of the perturbed eigenfunctions as well as the unitarity 
of the time- evolution operator |12j . 
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